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A NEW CONCEPT OF FLOW IN ROUGH CONDUITS 


Henry M. Morris 


SYNOPSIS 


A new and expanded concept of flow over rough pipe and channel surfaces, 
based particularly upon the effect of the longitudinal spacing of surface rough- 
ness elements and their associated vorticity streams, recognizes three basic 
types of rough conduit flow. These basic flow types have been denoted as 
isolated-roughness flow, wake-interference flow, and quasi-smooth or skim- 
ming flow. 

The friction factor equations derived for these three types of flow appear 
to be adequately substantiated by data obtained by many different experimenters, 
on all types of roughnesses, in both pipes and open channels. They can be 
extended to surfaces of variable roughness, so that reasonably accurate pre- 
dictions of friction factor can be made for all types of conduit surfaces. The 
concept is physically reasonable and is relatively simple analytically, so that 
it appears to bear promise as a more useful and rational method for the hy- 
draulic design of conduits than methods in current use. 


Introduction 


As a contribution toward the practical solution of the important problem of 

the hydraulics of flow in rough conduits, a general concept of such flow will 

here be offered, which is intended to make possible the prediction of friction 

factors for any type of conduit surface. 2 
The Darcy equation for head loss is as follows: 


f D2 (1) 
where H, represents the feet of head loss in a length L of the conduit, f is the 


f 
Darcy friction factor, D is the equivalent diameter of the conduit, and V is the 
average velocity of flow across any section. The friction factor is a dimen- 
sionless coefficient, which can be shown by dimensional analysis to be a func- 
tion of the Reynolds number of the flow and of several dimensionless geomet- 
ric ratios. Thus: 


f = $(R, h/D, A/D, s/D, - - -) (2) 


where R is the Reynolds number DV/?, ¥Y is the kinematic viscosity of the 
fluid, h is the radial height of the wall roughness projections, A is the spacing 
of the wall roughness elements in the direction of flow, and s is the clear 
peripheral spacing between elements. If other dimensions are necessary to 
describe the roughness elements, then as many more dimensionless geomet- 
ric terms should be included in the function d . 

The velocity distribution curve in pipes with rough walls has been found 
experimentally, especially by Nikuradse (11), to fit the following: 


= a+b log.) (3) 
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where a is a dimensionless parameter depending on the type of roughness, 
to is the pipe radius, and b is 5.75. The constant 5.75 is actually —— ‘ 


where k is the von Karman universal turbulence constant, which has been 
shown to be applicable to the interior turbulence of smooth and rough-walled 
conduits alike, and to have an approximate experimental value of 0.40. The 
velocity traverses obtained by Nikuradse in his sand-coated pipes are usually 


described by the following equation: 


= 8.5+5.75 (log,, 2) (4) 
where e is the diameter of uniform sand-grains coated on the wall. The cor- 


responding equation for friction factor is: 


1 2 1.74 (5) 
vi ° 


This equation obviously implies f to be independent of R and to be a function 


solely of the relative roughness — . Thus, this type of flow, commonly known 
To 


as rough-pipe flow, follows the quadratic resistance law, with the head loss 
proportional to the square of the velocity. 

However, Nikuradse also showed that rough pipes normally will experience 
a brief regime of smooth-pipe flow at low Reynolds numbers before the rough- 
ness projections begin to have a material effect on the turbulence. A transi- 
tion flow regime is then established between smooth and rough-pipe flow, in 
which the friction factor first decreases, then increases and finally becomes 
constant as R increases. 

It has often been assumed that all types of conduit surfaces would also give 
a constant friction factor at sufficiently high values of R. However, most 
commercial pipe tests have yielded descending f-R curves, with some tests 
giving horizontal characteristics and very infrequently rising characteristics. 

In recent years, the practice of using the Nikuradse equivalent sand-grain 
diameter, K., as a measure of surface roughness, has been increasing among 


engineers. This term is defined as the sand-grain diameter which would give 
the same value of friction factor as the actual pipe, obtained by solving eq. (5) 
for K_ =e. 

s 


Since, however, most f-R curves for commercial pipes have been found to 
have descending characteristics throughout the experimental range, the for- 
mula of Colebrook and White (2) is often used to determine the value of K. ; 


or at least to define the transition regime for the surface being tested. “This 
formula is: 


= 2log’o + 1.74-2log,, (6) 


ve «10K. RV 
Tr 


fe) 
/K, 
This equation gives a curve which is asymptotic to both the Nikuradse smooth 
and rough-pipe curves. If implies that f decreases as R increases. The 


2 log To 


parameter yr * 10 K, has been called by Prandtl the generalized 
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resistance function and the parameter fie can be shown to be a wall 
Ss 


Reynolds number equal to v°K, . V32. 

The use of the Nikuradse and Colebrook-White equations, together with the 
design methods based thereon, undoubtedly represents a substantial improve- 
ment in previous design methodology. However, there exist several impor- 
tant deficiencies in this method which it is hoped the concept of the present 
paper will overcome. 

Probably the most important such defect lies in the failure of the Colebrook- 
White equation to apply to conduit surfaces which produce either horizontal or 
rising f-R curves. The tests of Harris (5) on surfaces of “random” roughness 
yielded f-R curves which were essentially horizontal after leaving the smooth- 
pipe curve. Tests at the St. Anthony Falls Hydraulic Laboratory on commer- 
cial corrugated metal culvert pipes (16) revealed a definitely rising charac- 
teristic for the f-R relation. 

Furthermore, even for pipes which give the usual falling characteristic, 
experimental values of K. do not always turn out to be a function of only the 


type of surface. Tests on concrete pipes at the St. Anthony Falls Hydraulic 
Laboratory, for example, (14), (15), have implied K. as computed from the 


Colebrook-White equation to vary with both R and To 
Finally, it is apparent that K. has no definite physical significance as far 


as the roughness elements on the surface are concerned, with the exception 
of uniform sand-coated surfaces. However, its use has often seemed to imply 
that it is being associated with the height of projection h of the elements. 
There is no reason, of course, why h and K, should be equal or even be cor- 


related with each other, since other roughness dimensions will obviously exert 
a material effect on the friction factor. 

It is thus desirable, if possible, to develop a new concept of rough conduit 
flow which will utilize a physically realistic measure of surface roughness, 
and which will be able to explain and predict the different forms of f-R curves, 
and finally which can predict friction factors for any surface of known rough- 
ness dimensions and for any value of Reynolds number. 


The Importance of Form Drag and Roughness Spacing 


The chief source of friction loss in a fluid flowing over a rough surface is 
evidently the generation, spreading, and subsequent dissipation of vortices 
from the wake and separation zone behind each roughness element. Each ele- 
ment is thus a source of vorticity, and the longitudinal frequency with which 
such sources of vorticity occur is therefore of greatest significance in con- 
trolling the turbulence structure and energy dissipation phenomena. 

This fact strongly suggests that the longitudinal spacing ) of the roughness 
elements is the roughness dimension that is of greatest importance in rough 
conduit flow and would thus offer the most likely basis of correlation of fric- 
tion loss data in such conduits. 

With this in mind, two parameters involving the roughness spacing ) will 
be defined for later reference. The roughness index is defined as d/h, the 
ratio of the roughness element spacing to the radial height of projection. The 
relative roughness spacing is defined as To/rh? the ratio of radius to rough- 


ness element spacing. This term is similar to, but should not be confused 
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with the ratio usually defined as the relative roughness, Tr, /h. The radius 


lg, aS well as other radial dimensions y , will be measured with reference to 

a datum through the crests of the roughness elements. This is believed to be i 
the best choice of datum in view of the fact that the wake and vortex- generating 
zone at each roughn?ss element originates at or near the crest of the element. j 


Basic Types of Flow in Rough Conduits 


Recognizing, then, the fundamental significance of the roughness element 

Spacing in the expenditure of flow energy, it next becomes apparent that three 

different types of flow phenomena may occur near the wall; the type which 

occurs in a given situation will depend primarily apon the value of A for the 

given roughness. 
If the wall roughness elements are far apart, separated by long reaches of 

relatively smooth wall surface, then the individual elements will act as isolated | 


bodies on which are exerted drag forces by the flowing fluid. The wake and 
vortex-generating zone at each element is completely developed and dissipated 
before the next element is reached. The apparent friction factor would there- 
fore result from the form drag on the roughness elements in addition to the 
friction drag on the wall surface between elements. The form drag on any 

one element would depend largely on the height of projection h of the element. 
The total form drag in a given length of a conduit would depend upon A . Thus, 
the roughness index A/n, could be expected to be a significant correlating 
parameter, for this type of flow, which canbe designatedas isolated- roughness 
flow. 

The second type of flow results when the roughness elements are placed 
sufficiently close together so that the zone of separation and vortex generation 
and dissipation associated with each element is not completely developed be- 
fore the next element is encountered. This type of flow may be called wake- 
interference flow, involving as it does the interference by the wake of each 
element with the vortex spreading and dissipation connected with the preceding 
element. The friction drag of the wall itself is obliterated, and the entire re- 
gion near the wall is the scene of intense and complex vorticity and turbulent 
mixing. The height h of the element is relatively unimportant in this type of 
flow, but again the spacing A is obviously of major importance. Also the pipe 
size will control in part the radial extent of this wall region of abnormal tur- 
bulence. For wake-interference flow, therefore, the relative roughness spac- 


ing *o/X will be an important correlating parameter. 

Still another type of flow would result when the roughness elements are so 
close together that the flow essentially skims the crests of the elements. In 
the grooves between the elements will be regions of dead water containing 
stable vortices. The bulk flow will be over a pseudo- wall composed of the 
roughness crests and the upper limbs of the groove vortices. Large rough- 
ness projections will be absent from this pseudo-wall, and the flow will thus 
be similar to smooth-conduit flow. This type of flow might be called either 
quasi-smooth flow or skimming flow. Much of the energy loss in such flow 
can probably be attributed to the maintenance of the groove vertices, the size 
of which is controlled by the width or the depth of groove, whichever is small- 
er. The number of such vortices depends on A. The roughness index, ex- 
pressed either as A/h or rj, where j is the groove width, can thus be 
expected to be a significant parameter influencing the apparent friction factor 
in skimming flow. 
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Theory of Wake-Interference Flow 


The wake-interference flow phenomena will be considered first, primarily 
because of the fact that the well-known Nikuradse sand-coated pipes seem to 
have produced this type of flow. 

For the purpose of developing general equations descriptive of wake- inter- 
ference flow, consider a circular pipe section, at a distance sufficiently re- 
moved from the entrance or other irregularity to permit full development of 
the turbulent boundary layer. At each roughness element will be formed a 
wake, from which vortices will be shed towards the pipe center. The forward 
velocity of flow will wash the vortices downstream, where they will commingle 
with the vortices from the next roughness wake. Thus a relatively distinct 
wall zone of flow will result, in which will prevail a high comparstive intensi- 
ty of turbulent mixing. 

In the central region of the pipe, normal turbulence, characterized by the 
von Karman similarity constant k , will prevail. The boundary between the 
central and wall regions will be at a distance from the wall which is probably 
primarily dependent on the spacing of the elements, since it is this parameter 
which controls the frequency of vorticity sources and the spread angle of vor- 
tices toward the central region. Accordingly, the distance of the boundary be- 
tween the wall and central zones will be taken as cA, where c can probably 
be expected to vary somewhat with the wall Reynolds number and with the 
relative roughness spacing. 

The general dimensionless equation of velocity distribution in the central 
region is then: 


v 1 a 1 


where v is the velocity at distance y from the roughness crests, and v* is 


the friction velocity. The dimensionless number a is the value of A at the 
boundary between the two regions at y =cA, andthe number A is the value 
of Ee at the distance y =A. Note that A is also (a - i log .c). 


Because of the higher relative intensity of turbulent mixing in the wall 
region, the velocity distribution would be expected to be more nearly uniform 
in this region than would correspond to the von Karman constant k. The 
equation of velocity distributions in this wall zone is not known of course, but 
might reasonably be assumed to be of the same general form as eq. (7) though 
with a different turbulence parameter replacing k. 

Accordingly, the following velocity distribution equation is postulated for 
the wall zone: 


~ a + Wlog, a Ay log. (8) 


Again, a is the value of A at ¥ = cA, the break in the velocit~ distribution 
and AY is the value of at y=A. The W corresponds 
"E in the normal velocity distribution equation (7). If = is plotted 

saalian log, } for a given pipe section, then % is the slope of the curve 


for the PPE a region and Wis the slope in the wall region. In general the 
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value of %, must equal or exceed the value of W, because of the abnormal 


turbulence in the wall zone and therefore the more nearly uniform velocity 
distribution (when plotted on a semi-logarithmic basis as indicated). 

The parameter W is evidently not a constant since it depends upon the 
wake and vortex-spreading phenomena in the wall region, which depend in 
turn upon the form and dimensions of the roughness elements (especially 
their longitudinal spacing) and upon the fluid velocity and viscosity. Fora 
given wall roughness, it is probable that Y will prove to be a function of the 


wall Reynolds number, defined as RY and of the cross-sectional form of 


r 
0/ 
the individual roughness elements, and that, for a given form of element, YW 


will approach 24 as the wall Reynolds number, Ry increases. It is as- 


sumed that W does not vary significantly with y. This assumption ignores 
at least a small transition zone, both near the wall and near the central core 
zone, but the use of a constant, average, value of W throughout the wall zone 
will materially simplify the ensuing analysis and will not significantly affect 
the results. 

It is next desirable to obtain an expression for the average velocity of flow 
through the entire cross-section. This can be done by integrating the elements 
of discharge across the section, then dividing by the total area. Thus: 


| an J vrar 


* 
ok Vv Av Av 


| (r,-y) (a+ log, dy + | 


(r,-y) (ast log.y) dy} (9) 
cA cA 


Omitting the intermediate steps, this expression ultimately reduces to: 


Vv 


1 1 


2 1 r 

+= log. (10) 
k 
where V is the average velocity of flow over the cross-section. Replacing 
x by its equivalent, {’ the following results: 

r 


This equation may be described as the general resistance equation for axial 
wake interference flow. 

A similar equation may be derived for two-dimensional wake-interference 
flow in an analagous manner, the final equation being as follows: 


ve 


where b is the distance from the roughness crests to the channel center. 
The left-hand member of equations (11) and (12) will be referred to as the 

resistance function. The first term in the right-hand member is the lower 

limiting value of the resistance function. The second term in the right-hand 


member is the corrective term, which approaches and eventually reaches 


(12) 
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zero at high values of the wall Reynolds number, since W approaches vA and 


¢ decreases as increases. 


The term k is the universal turbulence constant of yon Karman, usually 
taken to be about 0.40, although there is considerable evidence that a more 
generally accurate value would be about 0.37 or 0.38. If 0.40 is used, the 


resistance function simplifies to 2 


The term A is asserted here to be another universal constant, related to 
the wall zone turbulence. This has not been recognized heretofore and will 
be seen to have very significant consequences. It is obvious that A is con- 
stant for a given conduit and roughness type, since it is merely the value of 
Fe in the core velocity distribution at y = A, and since te at any constant 
value of y is independent of R and Ry; The magnitude of A may easily be 
scaled off the dimensionless velocity distribution plot at the ordinate y=). 

Consider now the situation for some other arbitrary roughness, with the 
same value of A as the first roughness but otherwise with different dimensions 
and shape. Now assume that the flow velocity is adjusted to the point where 
the phenomena of separation on the roughness elements are mechanically 
similar to the corresponding phenomena for flow at some given value of R or 
Ry on the first roughness. In general this would require a different value 


of Ry than used for the first roughness, but it is apparent that such a similar 


condition could be produced. The location of the separation point on the 
roughness elements can be adjusted by adjusting R. The turbulence phenomena 
in the wall zone (beyond the roughness crests, of course) depend primarily and 
directly on the character of the wakes rather than the form of the elements 
producing the wakes. 

Assuming then that the wake and wall turbulence phenomena at some value 
of Ry, for the second roughness are made mechanically similar to those for 


the first roughness at some other value of Ry the =. distributions for the 
two situations must be similar in the wall zone. The values of W, c, and AY, 
‘would therefore be equal in the two situations, and this would also mean that 
the Se distributions in the core would be the same. In other words, the 
plots of log = vs. ~. for the two situations would coincide throughout. 


Therefore, the value of A can be made the same for the two roughnesses. 
But since, as we have seen, A is constant for any one roughness, it must 
therefore have the same constant value for both roughnesses. A similar 
argument could be developed for roughnesses of different A - values. 

The inference, therefore, is that the term A is independent of surface 


roughness. Its magnitude may be determined by scaling Es from any core 


velocity distribution curve, at a distance y = A from the roughness crests. 
This in turn means that the limiting value of the resistance function, 


i. (A - 2) is also a universal constant, independent of surface roughness. 


Its magnitude, as well as the magnitude of A_ will be determined shortly from 
the available experimental data. 
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First, however, the implications of the additive term should be noted. In 
the regime in which the corrective term has not yet vanished, its magnitude 
continuously decreases as the wall Reynolds number increases. An analysis 
of equations (11) and (14) will show that this in turn means that for wake- 
interference flow in any given pipe the friction factor must increase with in- 
creasing Reynolds number, becoming and remaining constant when the cor- 
rective term vanishes. 

The in the case of axial conveniently be sim- 


plified to -w) A) since the term > would almost always be 
To 
negligible compared to the term 2 (SA). For a given value of the wall 


Reynolds n«mber, the - yw) will depend only on the form of 
roughness element. The term aa is simply a ratio of two radial dimen- 


sions, the distances from the roughness crests to the break in the velocity 
distribution curves and to the pipe center. It is reasonable to think that these 
dimensions would be geometrically similar for a given roughness form, and 
thus that their ratio would be independent of the relative roughness spacing 


x° 

Consequently, the corrective term may be considered as primarily a func- 
tion of the wall Reynolds number Ry» and the cross-sectional form of the 


roughness elements. This means that, if the resistance function 4 
R VE 


r 
o/ 

a single curve should result for each form of roughness element. Further- 
more, the curves for all forms of roughness will approach and reach the same 


), 


r 
OF - 2 log, 0 > ) is plotted against the wall Reynolds number, ( 


limiting value of the resistance function = (A - 3). The value of Ry at 


which the corrective term vanishes would be expected to depend on the form 
of the roughness elements. 


Velocity Distribution Data in Wake-Interference Flow 


There are available several sets of experimental velocity distribution data 
on conduits of regular surface roughness pattern. However, most of these test 
conduits were of rather small size. 

The only applicable velocity distribution data on large pipes of which the 
writer has knowledge are a series of measurements made at the St. Anthony 
Falls Hydraulic Laboratory on 24 inch and 36 inch diameter corrugated metal 
pipes (9). The two pipes were each 193 feet long. In all, some fifteen velocity 
distribution curves were obtained at sections at which the turbulent boundary 
layer was fully a A few typical curves are shown on Fig. 1. 


When plotted with te 


all indicated a definite break between the wall and core velocity distribution 
curves. From these plots, the values of k, W, and A were measured and 
plotted with the wall Reynolds number as common abscissa. 

The value of A was found to vary only slightly from its average value of 
8.70, the range of values being from 8.4 to 9.3, and the standard deviation 
0.24. In view of all the possible causes of data dispersion, this appears to 
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as abscissa and 108149 1 as ordinate, these curves 


| 


demonstrate satisfactorily the essential constancy of the number A for cor- 
rugated pipe at a value close to 8.70. 

Similarly, the values of k were found to be fairly constant, with an aver- 
age value of 0.36, and standard deviation of 0.03. The scatter of values of 
both A and k about their means was essentially random. 

However, it was evident that the parameter W increased with increasing 
values of R,- The scatter of the data was such that any curvilinear relation- 


ship could not be detected. The straight-line regression equation for this 


relation was: 
R VE 
155,000 (13) 
It is probable that the assumed straight line would curve over at high val- 
ues to become tangent to the line W= : - The minimum value of R,, at 
which this would occur could be computed by extrapolation of eq. (13) to the 


w = 1.00 + 


value of W = ; e. * 2.77. This computation yields the very high value of 


274,000 for the smallest value of _ at which the limiting resistance function 


might be reached. Since the maximum experimental value attained for the 
latter was 67,000, it is obvious that very high velocities would have to be 
reached before a constant condition could be established. 

This means also that the friction factor must continue to increase as the 
Reynolds number increases, to very high values of the latter, for corrugated 
surfaces. It will be seen shortly that experimental results of friction tests on 
these pipes confirm this implication. 

Turning now to consideration of other pertinent velocity distribution data 
that have been obtained by various experimenters, results from the work of 
Nikuradse (11), Schlichting (13), and Vanoni (19), will be briefly examined. 

The tests of Nikuradse on pipes roughened by densely-packed uniform 
send grains have hitherto been interpreted in terms of the pope roughness, 
However, it is that the relative roughness spacing was nu- 
merically equal to: — in his pipes, so that either parameter could be used 


without affecting fe of his numerical results. The Nikuradse rough- pipe 
velocity distribution equation can therefore be written: 


Vv 1 
= 8:5 + log, (14) 


A 


The value of A for surfaces densely packed with uniform sand grains is there- 
fore about 8.5. Thus, two pipe surfaces as radically different in character as 
corrugated metal and sand give practically identical results for the general- 
ized velocity distribution equations. This is strong confirmation of the signifi- 
cance of the roughness spacing in the flow structure. No such correlation as 
this is at all possible on the basis of the roughness height. 

The tests of Schlichting included many different types of roughness elements, 
instailed on test plates which were arranged to form one side of a rectangular 
conduit in which two-dimensional flow was obtained near the transverse center 
of the plate. Most of his test plates produced isolated- roughness flow, and 
will be mentioned later, in the discussion of that type of flow. Three plates, 
however, evidently produced wake-interference flow. 
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One of these was a plate coated with dense uniform-grain sand, 0.135 cm. 
in diameter. Another was a plate to which were attached spheres 0.41 cm. in 
diameter and with a spacing of 0.60 cm. in both directions. The third plate 
had angles extending transversely throughout the plate width, 0.30 cm. high, 
and spaced at 2.00 cm. intervals. 

Schlichting published velocity distribution curves for all his roughnesses. 
Though the tests were all made at high values of Ry making the breakpoint 


in the velocity distribution curves difficult to determine closely, the latter all 
revealed a definite decrease in slope in the wall zone. 

After correcting Schlichting’s data and curves to a datum at the roughness 
crests, then determining A as the value of z. at y =A, it was found that 
the sand roughness plate and sphere roughness plate gave 7.8 for A; the plate 
with long angles gave 8.2 for A. These values are not very different from the 
A-values obtained in the Nikuradse and St. Anthony Falls pipes. They are 
consistently lower, however, which probably reflects the two-dimensional 
structure of the flow. 

The velocity distribution data of Vanoni were obtained in a rectangular 
open channel, with a bottom surface composed of densely packed sand, 0.88 
mm. in diameter. The published velocity profiles were taken at the center 
line, where the flow was two-dimensional. Those which are pertinent here 
consist of the three runs made with clear water (most of his data involved 
silt- laden flows). 

The boundary between the center and wall zones was quite distinct, occur- 
ring at a distance of about y = 10A. This large value of c (c was of the 
order of magnitude of unity for the St. Anthony Falls pipes and the Schlichting 
conduits) appear merely to confirm the deduction that c should be large for 
small values of A/r,. 


The A - values for Vanoni’s three curves are found to be approximately 
8.0, 6.8, and 7.7, with an average of 7.5. This is very near to the values of 
A in two-dimensional flow found for Schlichting’s three test plates, the slight 
difference perhaps being attributable to free surface effects. 

The available velocity distribution data, though limited, thus seem to sup- 
port the previous theoretical inference that the term A should be a universal 
constant, independent of the type of wall surface. A small difference due to 
cross-sectional form seems to exist, however. For axial flow, A is about 8.7 
and for two-dimensional flow about 7.8. 

Numerical values for A andk can now be inserted in the generalized 
resistance equations (11) and (12). If k is assumed to be 0.40, for simpli- 
city, and A is taken as 8.7, and 7.8 for axial and two-dimensional flow, re- 
spectively, these equations become, respectively: 


- 2 = 1.75 + (<a) (2.5- v) (15) 


175 (£4) (16) 


4 b 


b 


It is very interesting, though probably fortuitous, that the limiting value of 
the resistance function in both axial and two-dimensional flow is 1.75. Since 
most other cross-sectional shapes could probably be referred to a circular 
section through use of the hydraulic radius (r, = 2R,,)) the number 1.75 
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_ 


} 

| 


appears to be a most remarkable universal constant, independent of both sur- 
face roughness and cross-sectional form. 


Experimental Verification of the Generalized Resistance Equation 
for Wake-Interference Flow 


It has been seen that such velocity distribution data as are available for the 
wake-interference type of flow support the general theoretical relations devel- 
oped previously. However, these data are limited, and much stronger confir- 
mation of the theory can be obtained from the available friction loss data in 
conduits with surface patterns of roughness producing this type of flow. 

Eqs. (15) and (16) indicate that, for wake-interference flow, the resistance 
function decreases from some maximum value to a minimum value of 1.75 as 
Ry increases. This transition would be expected to have a form depending 


upon the form of roughness elements, but the limiting value of 1.75 would be 
applicable to all types of roughness elements. The maximum possible value 
of the resistance function could be estimated by assuming a W of zero, anda 


value of 4 of unity. These values would give 3.52 as the maximum possible 


value of the resistance function (2.64 for two-dimensional flow). The maximum 
value for any particular type of roughness may be less than this, but it cannot 
be greater. The equation further implies that the friction factor will increase 


r, 
as the Reynolds number increases, for any given value of = . Also, for a 
given Reynolds number, the friction factor increases as _2 decreases. 


The foregoing description of the transitional zone and limits for wake- 
interference flow comprises a rigorous set of specifications which experi- 
mental data must fulfill. However, it will be seen that every feature of these 
specifications is supported by a great variety of data, obtained by many dif- 
ferent experimenters, working under many different conditions and with many 
types of conduits, including both pipes and open channels. 

These data can best be analyzed when plotted with the resistance function 
as ordinate and the wall Reynolds number as abscissa. A logarithmic scale 
is used for the latter, because of the large range of values involved. The 
ordinates are plotted increasing downward in order to emphasize the similar 
appearance of this type of plot to the more familiar f - R plotting. 

First the data available for corrugated strip roughnesses will be examined. 
Tests on several commercial corrugated metal pipes have been carried out 
at the St. Anthony Falls Hydraulic Laboratory (14), (16). The tests included 
18', 24", and 36" diameter circular pipes and 18", 24'', and 36' pipe arch 
sections, The corrugations had a spacing A of 2 2/3 inches and height h of 
1/2 inch. 

It can be shown that these data strikingly exhibit the phenomena associated 
with wake-interference flow. Fig. 2 shows the resistance function plotted 
against the wall Reynolds number for these tests. The strong conformity of 
the data from all pipes to one curve is evident. The resistance function is 
seen to decrease gradually from about 2.7 to about 1.9, over a range of Ry 


from about 7,000 to60,000. The limiting value of 1.75 was not reached, but the 
curve gives every indication of a tangential approach to this limit as a value 
of Ry of probably about 350,000. This is essentially the same inference 


deduced from the W - Ry curve for the corrugated pipe velocity distributions 
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discussed previously. The friction data for part-full flow in these pipes were 
also found to lie on the same curve although with somewhat more dispersion. 

There are a few other sets of experiments on corrugated pipe surfaces of 
importance here. Gibson (4) made tests on a 2-inch copper pipe, with sine 
wave corrugations of ) = 0.4 inch and h = 0.1 inch. Streeter (17) made tests 
on two corrugation type roughnesses, numbered IV and VI in his notation. 

The tabulations of data published by Gibson and Streeter had to be recomputed 
on the basis of a datum at the roughness crests, and this was also true of 
most of the other published data on other types of roughness. 

The combined data from the results at the St. Anthony Falls Laboratory 
and those of Gibson and Streeter are plotted on Fig. 3. They define almost 
the complete transition curve from smooth pipe flow to quadratic-law turbu- 
lent flow. The curve is valid for a wide range of values of A, h, and Ty for 


both circular and deformed sections and for both pipe and open-channel tran- 
quil flow. 

Other tests on corrugated metal pipes, such as those performed by Yarnell, 
Nagler, and Woodward (20), did not include measurements of water tempera- 
ture, making it impossible to compute accurate values of R or Ry: However, 


the range in resistance function was about the same as in the St. Anthony Falls 
Tests. 

Fig. 4 contains resistance function data based on experiments on surfaces 
containing uniform spot-type roughness elements. These data include the 
familiar Nikuradse sand surface plot, the Schlichting sand roughness and 
spherical roughness previously described, and two roughnesses tested in an 
open channel by Bazin (1) and discussed by Keulegan (8). One of these rough- 
nesses was composed of fine gravel, the other of coarse gravel. 

Fig. 5 shows data for sharp-edged strip roughnesses. These include data 
obtained by Streeter (17) on two strip roughnesses composed of threads cut 
in a 2-inch brass pipe. One of these was trapezoidal in cross-section with 
d = 0.0435 inches and h = 0.012 inches (Roughness V in his terminology). 

The other was rectangular, with A = 0.087 inches and h = 0.022 inches. Also 
included are data obtained more recently by Streeter and others at Illinois 
Institute of Technology (18) on two square threac roughnesses in a 4.5 inch 
aluminum pipe. Schlichting’s angle roughness plate, referred to previously, 
is included. Also a surface tested by Fage (3) composed of staggered, 
square-based pyramids, is applicable here, since the pyramid edges formed 
what was essentially a sharp-edged strip roughness pattern. 

Finally, the approximate curves for all the data plotted on Figs. 2 through 
5, have been brought together in one plot on Fig. 6. This assemblage of curves 
strikingly delineates the nature of wake-interference flow, confirming all the 
theoretical implications previously noted. Each major roughness type pro- 
duces a resistance function curve which departs from the smooth-pipe curve 
at a value of the resistance function equal to or less than 3.5, and then gradu- 
ally decreases to the minimum value of 1.75, remaining constant thereafter. 
The limiting value is reached soonest for spot roughnesses, next for sharp- 
edged strip roughnesses, and las! for rounded strip roughnesses. 

Thus, the theory of wake-interference flow as developed previously, appears 
to be adequately substantiated by the available data. Furthermore, the curves 
of Fig. 6, judiciously applied, may also be used to predict friction factors for 
this type of flow in conduits of almost any regular roughness pattern. 
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Isolated-Roughness Flow 


This type of flow is characterized by independent development and dissi- 
pation of form turbulence behind each element, and by friction drag turbu- 
lence along the wall between the elements. The over all pattern of flow will 
be intermediate between smooth- pipe flow and flow following the square law. 
The friction factor would decrease with increasing Reynolds numbers, and in 
some instances might approach or reach a constant value at high Reynolds 
numbers. 

As a reasonable approximation, the total drag force may te written as a 
summation of wall friction drag and roughness element form drag, as follows: 
Ap(A) = 7,L, + 1/20Cpv,? h (p - ns) (17) 
where Ap is the pressure drop in length L, A is the area of flow cross- 
section, Ty is the wall shearing stress, p is the wetted perimeter, P the 


fluid density, h the roughness element height, Cy the drag coefficient for the 


element, n the number of individual elements ina periphery, s the clear trans- 


verse spacing between elements, and Ve the flow velocity at the roughness 


crests. 
If f is the bulk friction factor and i. is the friction factor for the wall sur- 


face alone, this equation can easily be shown to reduce to 
4C 2 
D ns 

( as) (18) 

Assuming the wall surface to be essentially smooth, it has been shown by 
v 

Nikuradse that 
be introduced and equation (18) then becomes: 


81-2 Cp (- (19) 
P 


This equation indicates the nature of isolated roughness flow. The friction 
factor is obtained by multiplying the smooth pipe friction factor by a coeffi- 
cient greater than unity, the excess depending on the roughness index )_, 
h 
and the proportion of normal wall length intervening 


2 
= 16.8 i As an approximation, this relation may 


the drag coefficient Cy 


between roughness elements around the periphery to the total wetted perimeter, 
~ For strip roughnesses, ns = 0. 

As a general check on the validity of equation (19), reference will be made 
to the various test roughnesses of Schlichting (13). These roughnesses . 
included spheres (two sizes), spherical segments, cones, short angles and 
long angles, each tested for various values of _A_ ands. 

These plates, with four exceptions, produced isolated roughness flow. 
Three of them produced wake-interference flow and have already been dis- 
cussed. One produced skimming flow. The measured friction factors on the 
plates giving isolated roughness flow agree very satisfactorily with friction 
factors computed from eq. (19). 

Considering also this type of flow in open channels, reference may be 
made particularly to the data of Powell (12), who tested flows over various 
arrangements of 1/4 inch and 1/8 inch square battens on the bottom and sides 
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of a rectangular flume. Again a comparison of measured friction factors 
with values computed from eq. (19) indicates excellent agreement. It is in- 
teresting that this agreement is obtained for both tranquil and rapid (though 
not “ultra-rapid” ) flows. 

The detailed tabulations of these computed and observed values of “f” are 
not included herein for want of space, but are available elsewhere (9) for 
those interested. 

It is believed, therefore, that the concept of isolated roughness flow 
harrmonizes with the data of Schlichting and Powell, and that eq. (19) there- 
fore gives a good representation of this type of flow. 


Quasi-Smooth Flow 


For surfaces with very low roughness indices, the flow will more or less 
skip over the crests of the roughness elements. Between the elements will 
be regions of dead water, with stable vortices set up in these regions by the 
moving fluid above the elements. The combination of the roughness crests 
and the still water in the interstices serves as a pseudo-wall, probably some- 
what wavy in character but without large roughness projections. Consequently, 
this type of flow has been called “quasi-smooth” or “skimming” flow. 

Vortex generation at the hypothetical smooth wall is the predominant 
source of turbulence in the bulk flow. However, maintenance of the stationary 
eddies between the elements probably requires expenditure of more of the 
available energy. As an approximation to the energy balance involved, it may 
be reasonably assumed that the total energy consumed is equal to the energy 
that would be required to maintain flow in a smooth pipe plus that required to 
maintain the groove eddies. 

These eddies can be regarded as essentially rectilinear circular vortices, 
of diameter equal to j, the width of groove (or h, the depth, whichever is 
smaller) and of length equal to the conduit perimeter p. The velocity at the 
vortex perimeter is equal to a coefficient C., times the wall velocity of flow 


V,,- The angular velocity w is constant at any radius. 


~~ The energy of flow per unit time through any concentric cylindrical shell 
of the vortex is: 


1 
dE = (dm) v2 = 3 P p(ar)v 3 


dole 


a. 69 
2 


The total vortex energy is: 


E = 2 2 3 rdr = 8 
j 

The head expended in a length of pipe L due to these groove vortices is 
therefore: 

Vv 8 AVv A 4R 

where Avy is the weight of fluid flowing per unit time. Adding the head loss 
ina smooth- pipe, and simplifying, the resulting bulk friction factor is therefore: 
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(23) 
(j,h) 


In the roughness index term, hn OF > would be used depending on whether 


h or j were smaller. The coefficient c. would have to be somewhat small- 


er than unity, perhaps of the order of magnitude of + The ratio would 


also be of this general order of magnitude. 

In many cases the additive term would be practically constant, and the f-R 
curve would be roughly parallel to that for smooth pipes (apparently diverging 
somewhat if plotted on a logarithmic scale, of course). 


Vv 
On the other hand, if there is a slight increase in a with increasing R, 


as seems possible, the additive term will be increased materially because of 
the cubing operation required. Consequently, for many roughnesses, the 
friction-factor curve in skimming flow may be expected to diverge somewhat 
from the smooth pipe curve at high Reynolds numbers. 

The best available experimental data on skimming flow appear to be those 
obtained at Illinois Institute of Technology under the direction of Streeter (18). 
Friction loss and velocity distribution measurements were made on three 
geometrically similar square-thread roughnesses cut into a 4 1/2 inch dia. 
pipe. The pertinent dimensions are given on Fig. 7. 

Note that the roughness index was the same for all three roughnesses, 
although the relative roughness was made to vary considerably. Eq. (23) 
indicates that the friction factor should be independent of the relative rough- 
ness, but should depend on the roughness index. Fig. 7, showing the f-R 
curves obtained on these pipes, shows the curves to be essentially coincident, 
just as implied by eq. (23). The curves also follow the implied trend of more 
or less paralleling the smooth-pipe curve, but finally diverging from it and 
evidently becoming essentially horizontal. 

The marked difference between these curves and those for the two other 
I.I.T. roughnesses (roughnesses II and III, shown on Fig. (5), may be noted 
in passing. The roughness elements for the two latter were identical in 
every way with roughnesses 1. with the sole exception of the spacing A , 


which was large enough to develop wake-interference flow, with its typical 
rising friction factor characteristic. 
The velocity distribution data ene on in the I.I.T. Report for roughnesses 


Ip _ , and ‘ seemed to indicate that — fluctuated randomly from about 0.62 


he aan 0.67, with an average of about a 65. Insertion of this value in eq. (23), 
together with various corresponding values of f and fy then reveals that the 


coefficient Cy must have averaged about 0.45, with small variation. This is 


clearly of the right order of magnitude, which gives further support to eq. (23) 
as providing a satisfactory skimming-flow equation. 

The only one of Schlichting’s roughness plates to produce skimming flow 
seems to have been his Plate V, which contained tightly packed spheres 0.41 
cm. in diameter. This plate gave materially smaller friction factors than his 
Plate II, on which the same spheres were placed at a 0.60 cm. spacing, and 
which has previously been noted to have produced typical wake- interference 
flow. 
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In spite of the differences between strip and spot roughnesses, it is inter- 
esting to examine the latter in boy light of eq. (23) as an approximation. It 


can be found that a value of 5 — of 0.65 and for Cc, of 0.52 will satisfy 


Schlichting’s experimental aati throughout the test range. These values 
appear quite reasonable. 

Though more data are needed, it therefore appears that eq. (23) may prove 
satisfactory as at least a preliminary quide to the character of skimming 
flow. 


Joint Roughness 


Special cases of isolated roughness and skimming flow may be due to the 
introduction of small-scale vorticity at the pipe or flume joints, If the joint 
causes a peripheral protrusion into the flow, then the isolated roughness 
equation (19) applies. If the joint causes a separation between adjacent con- 
duit sections, then the skimming flow equation (23) would be applicable. Both 
conditions would result in a descending f-R curve, approximately parallel to 
the smooth-pipe curve. 

This phenomenon is illustrated in two sets of experiments conducted at the 
St. Anthony Falls Hydraulic Laboratory. One was a 4-inch dia. Lucite pipe, 
with joints at 4.37 ft. intervals (10), the other on three large concrete pipes, 
18, 24, and 36 inches in diameter, each with a joint spacing of 6 feet (15). The 
joints on the 18! pipe were considerably smoother and more snugly fitting 
than on the 24'' and 36! pipes. 

All of these pipes were smooth-walled. The concrete pipes had been manu- 
factured by the cast-and-vibrated process, which produces a glossy-smooth 
surface. The f-R curves for these pipes are shown on Fig. 8, in each case 
being slightly above and essentially parallel to the Nikuradse smooth pipe 
curve. 

The joint characteristics of the concrete pipes plainly appear in the f-R 
curves. The 18-inch pipe has the lowermost curve, indicating the best joints. 
The 24-inch and 36-inch curves were essentially identical and at an approxi- 
mately constant increment of friction factor above the 18-inch curve. If the 
turbulence were being determined by the wall surface, as per the Nikuradse- 
Colebrook equations, the 18-inch curve should lie highest. 


Discriminating Criteria for Flow Types 


Equations have now been developed for each of the three basic flow types, 
but there is still the need of adequate criteria for determining which type will 
prevail in a given situation. 

The boundary between wake-interference and isolated- roughness flow can 
be denoted by equating the resistance functions for the two types of flow. 


Thus: 
1 


A/h p 


If the criterion is applied for a Reynolds number of say 2,500,000, for 
which f. = 0.01, and if the additive term in the right-hand member is assumed 


r. 
- 2 log = 2.9802 


1 cr 


to be negligibiy small, the equation simplifies to: 
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67.2 

100 

2 
(2 log _ + 1,75) 


This expression would probably prove syfficiently accurate as a discriminating 
criterion in most cases. The left-hand member involves the roughness element 
characteristics, the right-hand member the relative roughness spacing. For 
a given type of roughness element, and a given radius, the equation can be 
solved by trial for the critical value of the spacing ) which will mark the 
boundary between isolated- roughness and wake-interference flow. 

For example, solution of eq. (25) for Schlichting’s long angle roughnesses 
will be found to yield a critical A of 1.85 cm. Schlichting’s experimental 
data indicated that wake-interference flow was obtained when this type of 
roughness had a A of 2.0 cm., but isolated-roughness flow at higher values 
ofA. 

The spherical roughnesses of Schlichting gave wake-interference flow for 
a spacing of 0.60 cm., and isolated roughness flow at higher values. The 
equation indicates the limiting value to be 0.70 cm. 

Space does not permit detailed discussion of other data. However, it could 
be shown that the trends implied by the other available data on various types 
of roughness elements all are in substantial accord with the implications of 
eq. (25), which seems, therefore, to serve fairly satisfactorily as a discrimi- 
nating criterion between isolated- roughness and wke-interference flow. It 
may also be noted that, for all parameters except A held fixed, the friction 
factor increases as decreases in isolated-roughness flow. However, when 
wake-interference flow is established, a further decrease in A_ results ina 
decrease in friction factor. 

The criterion for differentiating between wake-interference and quasi- 
smooth flow cannot be set up in a similar way to that above, since the resist- 
ance function for quasi-smooth flow does not approach that for wake-interfer- 
ence flow as the A - value approaches the critical value between the two flow 
types, but rather moves away from it. The evident implication is that there 
is a sudden, rather than gradual, change from skimming to wake- interference 
flow. This sudden change evidently occurs when the stable vortex in the 
groove gives way to the typical flow separation phenomenon of a continual 
succession of vortices growing and washing downstream. 

It appears likely that this change will occur when the groove width j 
becomes significantly larger than the depth h. The vortex occupies the en- 
tire width of the groove as long as h is ‘arger than j and therefore can 
remain stable. Butif j is much larger than h, the vortex will adhere to 
the upstream face of the groove and the stream will flow over and down the 
vortex against the opposite groove face; this is essentially a typical flow 
separation action, and the vortex must rapidly become unstable. 

This inference appears to be upheld by the available data on skimming 
flow. The three I.I.T. roughnesses which produced skimming flow all had 
j/h ratios very slightly greater than unity. The same was true of Schlichting’s 
spherical roughness Plate V, which gave skimming flow. 


Flow Over Surfaces of Variable Roughness 


The formulas developed thus far have assumed a regular distribution of 
roughness elements of identical dimensions. A rational analysis is much 
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h (25) 
ns 


more difficult when the several parameters are permitted to vary. However, 
the three basic types of flow must still exist, and will occur when the rough- 

ness patterns satisfy the appropriate criteria. In many cases the best way to 
estimate friction factors for varied-roughness patterns would probably be to 

use simple statistical averages of the parameters which vary. 

Not many data are available as yet with which to check this assumption. 
However, the plank roughnesses of Bazin (1) can be considered as a special 
type of roughness with a repeating pair of A - values. These planks were 
equally spaced, but were much wider than their depth. This means that a 
wake would be produced at both the upstream and downstream edges of each 
plank, so that the average spacing of wakes was one-half the plank spacing. 
Bazin’s widely-spaced plank roughness consisted of 2.7 cm. by 1 cm. planks, 
on 7.7 cm. centers. This roughness produced the rising f-R curve charac- 
teristic of wake-interference flow, and the limiting resistance function, as 
computed for the mean spacing (i.e., one-half of 7.7 cm.), can be shown to be 
approximately 1.78. Similar results can be shown to apply for certain of 
Johnson’s plank roughnesses (7), which evidently produced wake-interference 
flow. The inference, therefore, is that an average value of A can be used for 
non-uniformly spaced roughness elements in the wake- interference flow 
equations. 

Spot roughness elements have been shown usually to produce isolated 
roughness flow, except when they are packed quite closely. Thus, each ele- 
ment acts independently and the friction factor could presumably be deter- 
mined by summation of appropriate drag forces on all the individual elements. 
However, when the roughness elements are closely packed, wake-interference 
flow is developed. If the roughness elements are predominantly of one size, 
with only occasional larger elements, the wake-interference phenomena will 
be controlled by the smaller elements, with isolated-roughness drag forces 
being added by the larger elements. 

These features are illustrated by the data of Colebrook and White (2), who 
tested various combinations of two sizes of sand coated on a pipe surface, the 
small grains being densely packed and the large grains in various isolated- 
roughness patterns. The small grains alone produced an f-R curve with the 
typical dip-and-rise obtained by Nikuradse and which implies wake-interfer- 
ence flow. The large grains alone gave the typical isolated- roughness de- 
scending f-R curve. 

The f-R curves for the various combinations seem to demonstrate that the 
isolated- roughness effects of the large grains controlled at low Reynolds 
numbers until the smaller grains came into play. Then wake-interference 
flow was produced, the bulk phenomenon giving a friction factor approximately 
equal to the friction factor for the densely- packed small grains plus the 
isolated- roughness value for the large grains. 

It can be seen from the Colebrook-White curves that, although wake-inter- 
ference flow was evidently produced, the dip and rise was considerably flat- 
tened by virtue of the isolated-roughness effect of the large grains. It seems 
likely that this effect would be enhanced by having more than one size of 
large grain interspersed among the small grains. The inference then is that, 
for densely packed spot roughness elements with a random size distribution, 
wake- interference flow would be produced but that the dip-and-rise would be 
practically eliminated, and that the f-R curve would become and remain hori- 
zontal very soon after leaving the smooth-pipe curve. This is, in fact, 
exactly the result obtained by Harris (5) in tests on pipes coated with stain- 
less steel spray, the resulting surfaces of which were similar to surfaces of 
densely-packed sand grains of randomly-varying sizes. 
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Types of Flow Over Commercial Conduit Surfaces 


It appears probable that the principles that have been delineated for flow 
over artificial roughness patterns can be extended to give a rational picture 
of the flow and resistance phenomena experimentally obtained on commer- 
cial pipe and channel surfaces. By far the majority of such pipes exhibit a 
decreasing friction factor curve, gradually diverging from the smooth-pipe 
curve. 

In many cases, this is attributable to joint roughness, as in the concrete 
pipes tested at the St. Anthony Falls Laboratory discussed previously. How- 
ever, most commercial pipe surfaces probably produce typical isolated- 
roughness elements such as rivets, projecting pieces of aggregate, etc. 
Frictionfactors for such surfaces should be predictable by means of the 
isolated- roughness flow equation (19), or in some cases by the skimming- 
flow equation (23) when the main roughness consists of a joint depression. 

Some pipes with very rough surfaces may produce wake- interference flow. 
In such cases, the typical dip-and-rise of the friction factor and resistance 
function curves will depend upon the form and degree of uniformity of the 
roughness elements. Uniform, rounded, elements such as corrugations will 
result in a deep dip and a long, gradual rise. Closely packed random rough- 
ness elements, on the other hand, will produce fully turbulent wake-interfer- 
ence flow almost as soon as the elements begin to pierce the laminar film, 
so that the friction factor curve will become and remain horizontal soon after 
leaving the smooth-pipe curve. 

Concrete surfaces have, in different tests, been found to give dropping, 
horizontal, and even rising, f-R curves. This variation obviously results 
from the wide variation of textures on different concrete surfaces. 


Summary and Conclusions 


The concept of flow over rough conduit surfaces which has been presented 
in this paper is based on the assumption that loss of energy in turbulent flow 
over such surfaces is attributable largely to the formation of wakes (and 
therefore vorticity “mills”) behind each roughness element. The frequency 
of such vorticity sources in the direction of flow will determine, to a large 
extent, the character of the turbulence and energy dissipation phenomena in 
the flow. Therefore the longitudinal spacing of the roughness elements is the 
roughness dimension of paramount importance in rough-conduit flow. 

Under this concept, there must exist three basic types of flow over rough 
surfaces. These flow types, together with their corresponding equations for 
friction factor, are as follows: 


1. Isolated-Roughness Flow: Flow over roughness elements sufficiently 
separated from each other so that the vorticity generated at one ele- 
ment is substantially dissipated before it reaches the next element. 


67.2Cq 
(19) 


The symbols have been defined previously, and are also given in the 
Glossary. For this type of flow, the friction factor decreases with 
increasing Reynolds numbers, and is independent of the radius tr 
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2. Wake-Interference Flow: Flow in which adjacent wakes and vortex- 
deterioration zones interfere with each other, causing a wall zone of 
abnormal turbulence. 


f = 1 (15) 


1 ch 
[2 logig + 1-75 + yr (2.5 - W) 
Oo 


The parameter W approaches 2.5 and c approaches zero, so that the 
last term in the bracketed pela vanishes at some value of the 


wall Reynolds number R th 


the roughness elements. This mears that the resistance function 
1 
1.75, for all types of roughness elements. When this regime is 
reached, the friction factor is independent of Reynolds number, and 


r 
is inversely dependent on the relative roughness spacing = . Before 


the exact depending upon the form of 


r 
10g will approach and reach the same constant value 


the constant value is attained, the friction factor increases as the 
Reynolds number increases. 


Quasi-Smooth (Skimming) Flow: Flow which essentially skims the 
crests of the roughness elements, with stable vortices in the grooves 
between the elements. 


In this type of flow, the friction factor decreases with increasing 
Reynolds numbers, and is independent of the radius ro 


The equation for determining whether isolated-roughness flow or wake- 
interference flow prevails in a given situation is as follows: 
» 
h 67.2 (25) 


____ 


D To 2 
(2 1.75) 


If the left-hand member exceeds the right-hand member, isolated- roughness 
flow will occur; conversely, wake-interference flow will occur if the right- 
hand member exceeds the left-hand member. 

The boundary between wake-interference flow and skimming flow is deter- 
mined by the groove dimensions. When the groove width becomes significantly 
larger than its depth, then the skimming action will change to wake-interfer- 
ence phenomena. 

The concept and the above equations may also be extended to surfaces of 
variable roughness, usually by merely employing average values of the rough- 
ness dimensions which vary. If the surface has roughness elements of such 
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variety that more than one of the basic flow types would be produced simul- 
taneously, then the friction factors for each flow type can be added to give 
the apparent friction factor for the surface as a whole. 

With regard to typical commercial conduit surfaces, it is apparent that 
isolated- roughness flow occurs over most of them, including concrete, rivet- 
ed steel, wood-stave, cast-iron, brick, etc., including also conduits which 
are nearly smooth but have protrusions at the joints. 

Wake-interference flow is obtained on corrugated surfaces, sand or gravel- 
coated surfaces, and probably in some old cast-iron and concrete pipes. 
Densely- packed “random roughness” elements will produce wake-interference 
flow which is fully developed as soon as the “smooth” flow stage is passed. 

Skimming flow is obtained in pipes or channels which are nearly smooth 
but have depressions at the joints. 

In general, the expanded concept of rough conduit flow presented in this 
paper appears to be quite adequately substantiated by experimental data from 
many different sources. Since it is both comprehensive in scope and also 
physically reasonable, as well as relatively simple analytically, it is hoped 
that it may serve as a stimulus to research and to the development of more 
rational practical engineering design methods in this field. It is suggested 
that, if the general concept and methods of the paper are adopted by the pro- 
fession, very useful corollary studies could be directed to the development of 
devices and methods for making rapid physical measurements of surface 
roughness characteristics, and tc the preparation of design charts or tables 
based on the various equations presented in the paper. 
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Glossary 


} 


re in the core velocity distribution at 


1 
(a - 4 log .¢) = the value of 


y=A. 


(a-wW log .c) = the value of 3 in the wall velocity distribution at 


y=A. 


the value of ‘. at y = cA, the breakpoint between the core and wall 
v 


velocity distributions. 

one-half the distance between the parallel flow boundaries in two- 
dimensional flow. 

the drag coefficient on an isolated roughness element, approximately 
the same numerically as for a body of similar size and shape in 
unconfined flow. 
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a coefficient such that the distance y = cA locates the breakpoint 
between the core and wall velocity distributions. 
a wall coefficient such that a. gives the tangential velocity of the 


upper limb of the groove vortex in quasi-smooth flow. 
the pipe diameter or equivalent diameter, ft. 

the uniform sand diameter in the Niku1adse tests, feet. 
the Darcy friction factor. 

the smooth- conduit friction factor. 


the friction factor due to the wall surface alone, exclusive of effects 
of roughness element protrusions or depressions. 
the loss of head due to normal internal friction in the flow, feet. 


the radial height of the wall roughness elements, feet. 

the width of the roughness grooves in skimming flow, feet. 

the von Karman universal turbulence constant. 

the equivalent sand-grain size, used as a measure of surface rough- 
ness, feet. 

the length of the conduit reach, feet. 

the number of individual roughness elements in a typical wetted 
periphery of conduit. 


= the wetted perimeter, feet. 


the conduit Reynolds number, 
the hydraulic radius, feet. 

R_ Vf 
the wall Reynolds number, 


the distance from the flow axis to any point in the flow cross-section, 
feet. 

the pipe radius, measured to the crests of the roughness elements, 
feet. 


= the clear circumferential spacing between adjacent roughness ele- 


ments in the periphery, feet. 

the average velocity of flow in the cross-section, feet per second. 
the velocity of flow at any point in the cross-section, feet per second. 
the velocity at the roughness crests, feet per second. 


the tangential velocity of any element in the groove vortex in 
skimming flow, feet per second. 


= the shear velocity,\/—>* , feet per second. 


= the radial distance from the roughness crests to any point in the 
cross-section, feet. 

the specific weight of the fluid, pounds per cubic foot. 

the slope of the dimensionless velocity distribution plot in the wall 


zone in wake-interference flow, that is, 


the longitudinal spacing of roughness elements, feet. 

the kinematic viscosity of the fluid, square feet per second. 
the fluid density, slugs per cubic foot. 

the wall shearing stress, pounds per square foot. 
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(17) 
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= the portion of the wall shearing stress due only to friction drag on the 


wall surface, pounds per square foot. 


= the angular velocity of the groove vortex in skimming flow, radians 


per second. 


References 


Bazin, H.: “Recherches hydrauliques,” Mem. divers savants, Sci. Math. 
et Phys., Vol. 19, Paris, 1865. 

Colebrook, C. F., and White, C. M.: “Experiments with Fluid Friction 
in Roughened Pipes”, Proc. Royal Soc. of London, Vol. 161, 1937. 
Fage, A.: “Fluid Flow in Rough Pipes”, Tech. Report, Aeron, Research 

Comm., Repts. and Mem. 1585, 1934. 
Gibson, A. H.: “The Flow of Water in a Corrugated Pipe” The 


Philosophical Magazine, Vol. 50, 1925. 
Harris, C. W.: “An Engineering Concept of Flow in Pipes” Trans, 


Amer. Soc. Civil Engineers., Vol. 115, 1950. 
Harris, C. W.: “The Influence of Random Roughness on Flow in Pipes”, 


Univ. of Washington Enging. Expt. Station Bulletins, Bulletin 115, 
1949. 
Johnson, J. W.: “Rectangular Artificial Roughness in Open Channels”, 


Trans. Amer. Geophys. Union, Section of Hydrology, 1944. 

Keulegan, G. H.: “Laws of Turbulent Flow in Open Channels”, Journ. 
Research Natl. Bureau of Standards, RP 1151, 1938. 

Morris, Henry M.: ‘A New Concept of Flow in Rough Conduits”, A 
thesis submitted to the University of Minnesota (Minneapolis, Minn.), 
in partial fulfillment of the requirements for the age of Doctor of 
Philosophy, 1950. 

Morris, Henry M.: ‘Preliminary Flow Tests on a Model Culvert”, Proj. 
Rept. 7, St. Anthony Falls Hydraulic Laboratory, 1948. 

Nikuradse, J.: “Stromungsgesetze in rauhen Rohren”, Forschungscheft, 
No. 361, 1933. 

Powell, R. W.: “Flow in a Channel of Definite Roughness”, Trans. Amer. 


Soc. Civil Engrs., Vol. 111, 1946. 


Schlichting, H.: “Experimentelle Untersuchungen zum Rauhigkeitsproblem”, 


Ingenieur-Archiv, Band VII, Heft. I, 1936. 

Straub, L. G. and Morris, H. M.: “Hydraulic Data Comparison of Con- 
crete and Corrugated Metal Culvert Pipes”, Technical Paper No. 3, 
Series B, St. Anthony Falls Hydraulic Lab., 1950. 

Straub, L. G. and Morris, H. M.: “Hydraulic Tests on Concrete Culvert 
Pipes” Technical Paper No. 4, Series B, St. Anthony Falls Hydraulic 
Laboratory, 1950. 

Straub, L. G. and Morris, H. M.: “Hydraulic Tests on Corrugated Metal 
Culvert Pipes”, Technical Paper No. 5, Series B, St. Anthony Falls 
Hydraulic Laboratory, 1950. 


Streeter, V. L.: “Frictional Resistance on Artificially Roughened Pipes”, 


Trans. Amer. Soc. Civil Engrs., Vol. 101, 1936. 


Streeter, V. L. and Chu, H.: “Fluid Flow and Heat Transfer in Artificially 


Roughened Pipes”, Final Report, Proj. 4918, Armour Research 
Foundation, 1949. 
Vanoni, Vito A.: “Transportation of Suspended Sediment by Water”, 


Trans. Amer. Soc. Civil Engrs., Vol 111, 1946. 
Yarnell, D. L., Nagler, F. A., and Woodward, S. M.: “The Flow of Water 


through Culverts” Univ. of Iowa Studies in Engineering, Bulletin 1, 
1926. 
390-23 


(1) 

| 

| 

| 


++ 


‘gee 


= = 


Fig. 1. Typical Velocity Distribution Curves (Corrugated Metal Pipes) A = 2 2/3 
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Fig. 2. Resistance Function (Corrugated Metal Pipes) 
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Fig. 3. Resistance Function in Wake-Interference Flow (Corrugation Strip Roughnesses) 
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Fig. 6. Approximate Summary Curves for Wake-Interference Flow 
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